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PREFACE 

The Agriculture and Resources Inventory Surveys Through Aerospace Remote 
Sensing is a multiyear program of research, development, evaluation, and appli- 
cation of aerospace remote sensing for agricultural resources, which began in 
fiscal year 1980. This program is a cooperative effort of the U.S. Department 
of Agriculture, the National Aeronautics and Space Administration, the National 
Oceanic and Atmospheric Administration (U.S. Department of Commerce), the 
Agency for International Development (U.S. Department of State), and the 
U.S. Department of the Interior. 

The work which is the subject of this document was performed by the Earth 
Resources Applications Division, Space and Life Sciences Directorate, Lyndon B. 
Johnson Space Center, National Aeronautics and Space Administration and 
Lockheed Engineering and Management Services Company, Inc. The task performed 
by Lockheed Engineering and Management Services Company, Inc., were 
accomplished under Contract NAS 9-15800. A major portion of this task was 
performed at Texas A&M University in partial fulfillment of requirements for 
a doctorate degree in statistics under professor H. 0. Hartley (supported on 
contrjct NAS9-13894). 
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1, INTRODUCTION 


1,1 THE AGRICULTURE AND RESOURCES INVENTORY THROUGH AEROSPACE REMOTE SENSING 


mmim] 


1,1,1 DESCRIPTION AND OBJECTIVE OF THE AgRISTARS PROGRAM 

The AgRISTARS program officially started in fiscal year 1980 after a success- 
ful completion of the Large Area Crop Inventory Experiment (LACIE). The LACIE 
was a project to test the technology that would be used in a satellite-aided 
crop production estimating system utilizing remote sensing. Instruments 
onboard Earth-observing satellites (i,e,, Landsat) were used to obtain infor- 
mation about the Earth by scanning its surface from the orbiting satellite. 


Over the past two decades, many developments led to the technologies used in 
LACIE. Some notable developments were: 

(1) Multispectral scanners capable of scanning the Earth's surface and 
producing a quantitative rj<Jiometric map at visible, near- infrared, and 
thermal infrared wave-1 engtiis; (2) pattern recognition techniques permit- 
ting agriculture crops to be identified on the basis of differences in 
spectral response during the year; (3) high-speed digital computers; 

(4) the 1972 launch by the National Aeronautics and Space Administration 
(NASA) of the first of the Land Satellite (Landsat) series of polar- 
orbiting satellites, making it possible to monitor each point on the 
globe every 18 days with multispectral scanners; (5) the development of a 
global weather-reporting network by the World Meteorological Organization 
(WMO); (6) the development of models capable of relating weather to crop 
yields (refs. 1 and 2); and (7) the development of many custom-made 
statistical methodologies. 

The LACIE completion proved that Earth-observing satellites used to gather 
data on agricultural and other resources can provide foreign agricultural 
production information with accuracy and timeliness. 


The successes of LACIE have stimulated an effort to "develop and evaluate 
aerospace remote sensing technology for other major commodities and global 
crop regions" (ref. 1, pp. 13-15) through its new program named AgRISTARS. 
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"The specific objectives of AgRISTARS Include development, testing, and 
evaluation of procedures for adapting space remote sensing technology to 
Improve the U.S. Department of Agriculture (USOA) capabilities to provide 
early warning and timely assessment of changes In crop conditions; to pro- 
vide more objective snd reliable crop production forecasts; to assist In 
Inventory and assessWitfnt of land, water, and other renewable resources; 
and to develop a cost base to permit the USDA to assess feasibility of 
Integrating space remote sensing technology with existing data systems." 
(ref. 1, pp. 13-15) 


1.1.2 GENERAL APPROACH OF THE PROJECT 

The project develops technology for use In processing Landsat data primarily 
to Inventory the amount and geographic distribution of crop acreage available 
for harvest. It uses meteorological data In agricultural -meteorological 
models to forecast yield per harvested acre (production). 

"This general approach requires that each crop region be stratified with 
relatively homogeneous subregions and randomly selected samples of Landsat 
data be machine processed to Identify and measure the areal extend of the 
crop types of Interest. Yield models developed over subregions by the 
Yield Model Development Project In AgRISTARS will be exercised to forecast 
yield. The estimates of area and yield will be combined to estimate area, 
yield, and production at the regional levels." (ref. 2, p. 1.4) 

Symbolically, the production estimation approach Is presented as follows. The 
crop region 1s stratified Into L strata: 1, 2, •••, h, •••, L; and 

(CROP ACRPAGE>5t^3t„ « (CROP - (CROP PRODUCTIOH),t^,j„„; 

(CROP PR00UCTI0N)^^g,„„ - (CROP PRODUCTION)st^gt,„ 
or (CROP PRODUCnON)^eglo„ > (CROP ACREAGE)^ 5 g,(,„ x' (CROP YIEL0)^jg,„„. 

1.1.3 SCOPE OF THE PROJECT 

The project Involves many different crops such as wheat, barley, rice, corn, 
soybeans, cotton, sorghum, and sunflowers, over many selected regions within 
the United States, Canada, India, U.S.S.R., China, Australia, Argentina, and 
Brazil . 
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SAMPLE DESIGN USED IN THE PROJECT 


1.1.4 

1*1. 4.1 Description of the Basic Sample Design 

The objective of a sample design for the project Is to estimate the crop 
acreage for each of J crops In a particular region of a specific country. The 
basic sample design of the project Is the stratified simple random sample 
design. At least, this basic design Is used In a ’Group I' situation*. Let 
Aj be the crop acreage for crop J, J * 1, 2, •••, 0 In a designated crop 
region of a specific country. 

The crop region is divided Into L subregions called strata which are relatively 
homogeneous In terms of the crop density (hence, the crop acreage, yield, and 
prod/^ctlon). 

The sampling units are the acreages of the 5-by-6 nautical mile area. These 
units are called 'segments'. 

For h * 1, 2, •*•, L, let Nh be the total number of segments made up of the 
h^^ stratum and n^ be the number of segments randomly selected by the simple 
random sampling (S'^S) frcm the N^ segments oi' stratum h. 

Lc>i be the crop acreage for crop J In the stratum h of the crop region 
under consideration and Aj,j^ be the total crop acreage for crop j In the 
segment 1 randomly selected In the stratum h. 

It 1s suppofied to have the mean acreage per segment as 

\j* “ 


*The ‘Group r situation is one ln v/hich the stratum sample size nu is greater 
than 1, h = 1, 2, •••, L, " 
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but the SRS plan gives only n^^ segments, 1.e., n|^ crop acreages 
1 « 1, 2, nu. However, the SRS design gives an unbiased estimate of J.., 
which Is 

V 

which Is the mean acreage of crop j per segment In the nj^ segments sampled In 
stratum h. 


The stratified unbiased estimator of Aj Is 





k 


^h^h j 


The unbiased estimator Aj has variance 


(3) 


Where S^j Is the variance within the h^*^ stratum of the segment acreage of 
crop j; I.e., 


where 



(5) 


( 6 ) 


the mean acreage of crop j ner segment in the N|^ segments of stratum h. 


Then, the estimate of the production of crop j In the region will be 


A 

(Production of crop j) region " ^ 


(7) 
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where 


Yj is the estimate of the average yield of crop j in the region under 
consideration. 


1.1. 4. 2 Flow Chart of a General Approach of Estimation 

A general approach used in estimation procedures can be seen in tl.e flow chart 
shown in figure 1-1. 


1.1. 4. 3 Some Drawbacks of the Stratified Simple Random Sample Design 

One important drawback of the stratified simple random sample design is that 
it does not utilize the fact of the existence of a significant positive 
correlation between the crop segment acreages observed in the current and 
previous years in order to yield more accurate estimates. 


Another shortcoming is that the stratified simple random sample design cannot 
overcome the problem of misping data. The stratified estimate 








which is shown as equation 3, will be badly effected when, for any reason, one 
stratum acreage estimate cannot be obtained. 


Moreover, it is supposed to choose each n|^ segment, h = 1, 2, L, out of 
Nf, segments of the h^*^ stratum randomly according to the simple random 
sampling plan. A list of all segments in each stratum h has to be created 
first, then a random procedure will help to choose n^, of is known in 

advance and n|^ is predetermined. If, for some reason, the randomness of the 
selection of sample segments cannot be preserved, the stratified simple random 
sample design may lead to some tragic results. 

Furthermore, the implementation cost of the stratified simple random sample 
design can be reduced by using the "rotation sample design" which is 
introduced in this paper. 

/ 
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Acreage estfmatlon 




Yield estlrnatlon 



Figure 1-1.- 


Flow chart of a general estimation approach. 



1.2 THE ROTATION SAMPLE DESIGN 


The rotation sample design which is presented in this paper was developed to 

(1) improve the accuracy of the crop acreage estimation procedures and 

(2) reduce the cost of the implementation. 

The objects of this design are not restricted to any particular crops but, for 
the sake of practical illustration and testing, wheat has been selected as the 
crop on which attention will be focused in this paper. 

1.2.1 OBJECTIVE OF THE ROTATION SAMPLE DESIGN 

The objective of the development of the rotation sample designs for the esti- 
mation of wheat acreages is to provide multiyear estimates of wheat acreages 
which are more accurate than stratified simple random sample estimators. The 
underlying assumption is that the variation of the wheat acreages of a partic- 
ular sample unit (named parcel or segment) from year to year is usually less 
than the variation of wheat acreages between segments within a particular 
year. 


1.2.2 GENERAL APPROACH 

The wheat acreage estimation in the design uses the available estimates of 
"matched" segments, i.e., segments for which estimates are available for two 
or more years. 

Instead of following strictly the simple random sampling plan to obtain np, out 
of N^ segments of stratum h, h = 1, 2, •••, L, for each year of estimation, 
the rotation design will retain each year a fraction of the sample segments 
and replace the remaining sample segments with new segments in accordance with 
a rotation pattern of those defined in section 3.2.2. 

The estimation theory is based on an additive-mixed analysis-of-variance model 
with years as fixed effects and parcels as a variable factor. 


1.2.3 OBJECTIVE OF THE ROTATION SAMPLE DESIGN TASK 


The objective of this effort is to develop a rotation sample design, survey- 
many rotation patterns, and highlight those rotation patterns for which the 
wheat acreage estimates for the current year have a minimum variance-reduction 
ratio. 
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2. STATUS 


Because of the specialized nature of this development item, the literature is 
very limited. Even if the term 'rotation design' is familiar to many statis- 
ticians in the field of sample surveys, especially to those who are familiar 
with the Current Population Survey (CPS) of the Bureau of the Census (ref. 3), 
the methodology of this particular rotation design is different from that of 
the CPS rotation design. For instance, the objectives for the AgRISTARS/FCPF 
project are considerably different from those of the CPS. 

2.1 THE BUREAU OF THE CENSUS CPS 

2.1.1 DESCRIPTION 

Of concern in the CPS is the estimation of social characteristics of the U.S. 
population and the provision of monthly estimates for these characteristics 
(i,e,, employment, unemployment, income distribution, family characteristics, 
marital status, migration, education, etc.). The estimates are made from a 
sample of segments of households which are arranged in "rotation groups" 
according to the CPS rotation design. There are eight rotation groups (eight 
systematic subsamples) of segments in each sample. A given rotation group is 
interviewed over an 8-month period which is divided into two 4-month periods. 
The segments in the selected rotation group are interviewed during the four 
consecutive months of the first year, then omitted from the survey during the 
following eight months, then interviewed again for the same four calendar 
months of the next year, and, finally, dropped from the survey. 

2.1.2 OBJECTIVES 

The basic reason for rotating part of the sample each month is to invalidate 
the problem of uncooperation of respondents. 

The other objective in using the CPS rotation design is to provide estimates 
which take advantage of accumulated information from earlier samples as well 
as the information from the current sample, which results in smaller variances 
of estimates. 
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2.1.3 PROCEDURE 


The estimator of a characterlstiTC, for example y, the so-called composite 
estimator, is a weighted average of the following two estimator components, 

a. The first component is the regular ratio estimate, for example y^, based 
on the entire sample for the month under current consideration time, t. 

b. The second component is the addition of the composite estimate for the 
preceding month t - 1, for example y^.p and the estimate of the change in 
each item from the preceding month to the current month. The estimates 
used to compute the change are the regular ratio estimates y^ and 
y^_l ^ (i.e., the segments that are in the sample in both month t - 1 and 
month t). The second estimate is then 




,t-l 


y 


t-i,t 


(81 


Then, the two estimates are combined as a weighted average, with weights 
summing to one. 


In summary, the CPS rotation design composite estimator is 


yj = W(yJ.j + . (1 - W)y; 


(9) 


where 

y!^ = the composite estimator for month t, 

y^ = the regular ratio estimator based on the earlier sample for month t, 

y^ t-1 “ regular ratio estimate- for month t but is made only from the 
matched segments which are in month t, 

y^_j^ ^ = the regular ratio ^'^tlmate for month t-1, but is made only from the 
matched segments wh'ch are in month t-1. 
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When there is a positive correlation between y values in two consecutive 
months, the variances of the changes ' ■^t~l,t^ small, since 

Hence, this rotation design with its composite estimator does utilize effec- 
tively the information from the earlier and current samples and results in 
smaller variances due to consideration of the estimates of change. 

2.2 A ROTATION SAMPLE DESIGN FOR THE AgRISTARS FCPF PROJECT 

The effectiveness of the composite estimator in a rotation design depends upon 
the year-to-year correlation between crop acreages. Indeed, there is usually 
a strong positive correlation between the wheat acreages of segments observed 
in consecutive years, and this similarity seems to imply that the experience 
gathered in the CPS possibly can be transferred for use in the AgRISTAilS FCPF ’ 
project. Such a transfer, however, .is not possible because of the differences 
in the objectives, the differences in segment construction (the "segment" in 
the CPS is a geographic segment comprising a number of households), and the 
s,'istantial differences listed as follows. 

a. The AgRISTARS FCPF time series is yearly and extremely short (two or three 
years). The CPS time series is monthly. 

b. The FCPF rotation sample design to be developed will have to be suitable 
to the peculiar trait whereby a considerable number of matched segments 
(i.e., segments in the sample for two consecutive years) will be lost 
through cloud cover or other reasons. So, the FCPF rotation design, in 
this case, has to provide estimators which are capable of dealing with 
unbalanced segment patterns over a moderate number of years. 

Henceforce, the CPS composite estimator is not suitable fo^ use in the 
AgRISTARS FCPF project. The mixed analysis of variance (ANOVA) models and the 
associated estimators provided by proper rotation designs were selected and 
combined to create a n'Svi' sample design which is custom-made for the AgRISTARS 
FCPF project. 
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This sample design will deal suitably with the completely unbalanced matching 
patterns and will utilize the fact that the variation of the wheat acreages of 
a particular segment from year to year is usually less than the variation of 
the wheat acreages of different segments within a particular year. 


3. PROCEDURE 


Before the main features of the procedure are introduced, it is worthwhile to 
review some basic background information. 

3.1 TARGET POPULATION AND SAMPLING FRAME 

3.1.1 POPULATION AND ITS OBJECTS 

Agricultural acreages over many selected reg'^ons within the United States, 
Canada, India, U.S.S.R., China, Australia, Argentina, and Brazil is the target 
population. Wheat is the crop of interest in this paper. However, any other 
crop such as barley, corn, or soybeans can take the place of wheat provided 
the crop characteristics satisfy the requirements of the procedure for wheat. 
One of the requirements is that the variation of the wheat acreages for a 
particular segment from year to year is usually less than the variation of the 
wheat acreages of different segments within a particular year, 

3.1.2 SAMPLING UNITS AND SAMPLING FRAME 

Elementary units are 'pixels' of the earth's surface. Each pixel is a picture 
element of the imagery provided by the multi spectral scanner (MSS) used by 
Landsat. It corresponds to an area of about 1 acre on the earth's surface. 

Sampling units are areal parcels which are clusters of pixels. The sampling 
unit used for che current AgRISTARS FCPF project is made up of 6- by 3. Simile 
segments, the imagery of which comprises approximately 22 932 pixels. 

The sampling unit size may be smaller than the 6- by 3.5-mile area as long as 
it supports the requirement of the procedure mentioned in section 3.1.1 and 
helps to reduce errors caused b^ classification, labeling, and calculation of 
the proportion of wheat within each segment. To explain more fully, classifi- 
cation, in automated analysis of remotely sensed data, is the process of 
assigning data points to feature classifications by a testing process in which 
an automated electronic system scans a total image (data set) pixel by pixel 
and determines whether the spectral properties of each pixel correlate with 
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those of the subject being classified. Labeling Is the process of determining 
characters which Identify an Item of data, an area of memory, a record, or a 
file. 

Each crop (wheat) region Is stratified Into a number of subregions, called 
strata, which are homogeneous In terms of the crop (wheat) density. This 
means crop (wheat) densities vary little within a stratum but may vary 
considerably from stratum to stratum. There Is a known number of strata, for 
example H, within a crop (wheat) region under consideration. Each stratum Is 
comprised of a known number, for example Nj, for stratum h, h « 1, 2, •••, H, 
of sampling units, e.g., 6- by 3.5-m11e segments. 

The sampling frame Is constructed by first covering a crop-growing region of a 
country by a large grid of segments that are either 6-by-3.5 miles or smaller, 
then excluding those segments which appeared to show less than 5 percent agr1= 
culture. The actual sampling segments (sampling units) are chosen from this 
sampling frame. 

3.2 ROTATION SAMPLE DESIGN 

The sample design presented In this paper is named the rotation sample design. 

3.2.1 SAMPLING PLAN 

After a crop (wheat) region is stratified into L strata and suppose stratum h 
is comprised of Nj^ segments, a well-defined allocation procedure will deter- 
mine the n|^ value, number of segments in stratum h to be in the sample. The 
nj, number will be the same for each year of study, but the segments will be 
chosen by a specified rotation pattern. 

3.2.2 ALLOCATIONS BY ROTATION PATTERNS 

With the present one-year-only design, the number of segments nj, allocated to 
stratum h varies from 1 to a sizable number such as 30 depending on the 
country of Interest. (The number n^ Is a number selected from Nj^ segments 
through a simple random sampling plan.) In the United States, the number is 
small; it is one, two, or three when counties are treated as strata. 
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The rotation sampling plan of this research is confined in the situation in 
which nj^ ■ 1, 2, or 3 segments to be in the samole under consideration in each 
year. In the situation named Group I, the stratified simple random sample 
design is taken as the current sample design of the AgRIStyiRS FCHF project. 

The allocation will be in a rotation pattern, as shown in figures 3-1 to 3-7, 

In figures 3-1 through 3-4, the notation S is the number nj, of segments per 
stratum in each year and the notation r is the number of years elapsing before 
a segment returns to the sample. Thus, figure 3-1 presents S » 2 segments per 
stratum, and each segment returns to the sample after r ° 2 years. 

In figures 3-5, 3-6, and 3-7, the notation r » « is used to indicate that a 
segment will never return to the sample. The expression inside the C ] shows 
the type of pattern, for example [2 retained 1] means that of the S ^ 2 
segments, one is retained in the next year. 

3.2.3 • SPECIAL CASES 

One special case occurs when some strata have only one segment included in the 
sample. In this case, the following principles of collapsing can be applied 
to match one of the specified patterns, i.e., one of those shown in 
figures 3-1 through 3-7. 

One stratum having only one segment included in the sample is collapsed with 
the other stratum having one of the following: 

a. one sample segment 

b. two sample segments 

c. three sample segments 

The collapsibility will be based on the relation of the stratum itself to the 
other strata which are alike. 
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Figure 3-4.- (S » 3, r “ 3) Patterns. 
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Figure 3-6.- (S =» 3, r = “) Patterns [3 retained 2]. 
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Figure 3-7.- (S = 4, r = ») Patterns [4 retained 3]. 


3.3 SPECIFIC APPROACH 

The mixed analysis of variance model is assumed in this rotation sample design 
in order to obtain estimators from the completely unbalanced matching 
patterns. The patterns are likely to arise through losses of segments (due to 
cloud cover or other reasons), and they differ considerably from any balanced 
rotation designs. 

This paper presents only the basic, simple model which will apply to many 
different rotation patterns. For each pattern, the variance of the stratum 
crop acreage estimate for the current year is derived, and an optimal design, 
which gives a minimum variance of the estimate, is obtained empirically. 


3.4 RESOLUTION 

3.4.1 MIXED ANALYSIS OF VARIANCE MODEL 

3. 4. 1.1 Basic Model 

The estimate of the average wheat acreage per segment is obtained from an 
infinite analysis of variance model for A^g, the "observed" wheat acreage of 
segment s in stratum h for year t. The model is written as 

\s " “t *^s ®ts 

for t = 1, 2, ••• , T years and s = 1, 2, ••*, sj^ segments 
where 

= average true wheat acreage per segment in year tj values are fixed 
year constants, 

bg = true segment variables applicable to all years with the assumption 
b ~ N(o,a2j 

e ^5 = composite segment error variable of segment s in year t with the assump- 
®ts ~ error variable contains two components which 

are the deviations of the true wheat acreage of segment s of stratum h 
in year t from the additive formula + bg and the classification error 
in Atg. 
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3.4. 1.2 Estimation 


3. 4. 1.2.1 Approach of Estimation for Average Wheat Acreage Per Segment 

The development of the basic model (11) to obtain the estimates and its 

variances is based on the assumption of a known constant y = '^b/'^e 
be estimated from pilot data. Therefore, the Aitken weighted least-squares 

A 

estimation will be applicable and utilized to obtain the estimators and its 
variances Var|a^j. 

The model shown in equation (11) can be written in the following matrix form: 

a = Xa + Ub + le (12) 

rtj fsj rsj y f 


where X and U are the design matrices of equation (11) and are to be defined 
according to each rotation design in the coming sections, and a is the 
vector of the observed wheat acreages of segments s in years t. In the model , 
the following values are held. 


a = (^ 11 ? ^1* ***> » ^12’ ***> ’ ***> ’ ***’ ^sh^* 

a = (a^» 02 » ••», «j)' 

b = (b^, b2, •••, bg, ..., bg^)' ; b ~ 


I = the identity matrix 

e = (e^]^, 621, •**’ ®ts> ®12> ***’ ®t2’ **•> ®ts> '**’ ®Tsh^' ’ S 



Since (a) is a vector of fixed-year constants, equation (12) gives 

Var(a) = Var(Ub + le) 


= U Var(b)U' + Var(e) 

= UU'o^ + u\ 
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(14) 


or, 


where 


0 


Var(a) = Ho^ 



H 5 I + yUU' 


(15) 

(16) 


with Y being regarded as a known constant, H is supposed to be invertible. 
Being multiplied both sides by equation (12) gives 

Ub + e (17) 

rw rw f\j 


Put 



H 

W 

H 


“l/^x = X* 

_ y* 





then equation (17) can be written as 

a* = X*a + U*b + e* 


where 

Var(a*) = Var|H"^^\j 

= Var(a)H’^^^ 



(18) 


( 19 ) 
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and 


Var(U*b + e*) = Var 


(h- 


1/2 


Ub + H 


■ 1 / 2 , 


= Var(b)U'H“^^^ + Var(e)H”^^^ 


= 


("“e - '4)" 


-1/2 ^ i^-l/2^2u-l/2 




because equations (13) and (14) imply 


U Var(b)U' « UU'o^ 

u 2 ,2 

= Ha^ - IC3 

Hence, equations (20) and (21) imply 

Var(U*b + e*) = H“^^^(H - I)H"^^^Og + 

= H"^/^(H - I + I)H“^^^a^ 

T 2 

= I^e 

Equations (19) and (22) give the following conclusion. 

Var(a*) = Var(U*b + e*) 

So, the model shown in equation (18), 

a* = X*a + U*b + e* 

has the following property 

Var (a*) = Var(U*b + e*) = l<? 

rs4 0 


(20) 


( 21 ) 


( 22 ) 

(23) 

(24) 


(25) 


3-10 


It 1s well-known linear model with known formulae which are deduced as 

follows. The normal equation of model (24) is 
* 

X*'a* = (X*'X*)a (26) 


That means 


i .e 


• > 




X'H“^a = X'H’^Xa 


The solution is 


a = (X'H"^X)"^X'H“^a 


(27) 


(28) 


which is the best linear unbiased estimator (BLUE) of a by the Gauss-Markoff 
theorem and is also the maximum likelihood estimator. From pilot data, y can 
be estimated as y = which is substituted in equations (27) and (28). 

If a consistent estimator of y is employed, it can be shown (see ref. 4) that 
equation (28) is still consistent and the variance of the estimator, which 

A 

will be given below in equation (29), still applies asymptotically with y Y 
as an error of order for magnitude Var|ajvar(Y) . Equation (28) gives 

Var(a) = Var[(X'H"^X)"^'H"^a] 

= (X'H"^X)"^X'H"^Var(a)H*^X(X'H"^X)‘^ 

= (X'H”^X)"^(X'H'^HH‘^X)(X'H'^X)"^ag 
= (X'H"^X)"^cTg 


that is. 



( 29 ) 
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3. 4. 1.2. 2 Approach of Estimation for Stratum Wheat Acreage 

From equations (28) and (29), the estimate of average wheat acreage per 
segment within stratum h at the current year Indexed t = T and Its variance 
and Var|aj|^j are obtained. 

Then the estimate of stratum h wheat acreage at the year T will be 
with Its varlanes formula 

Var(i,h) = <var(;rt,) 

where Var^o^j^j Is deduced from equation (29). 

3. 4. 1.2. 3 Flow Chart of a General Approach of Estimation 

As a general summary, the flow chart shown as figure 3-8 presents a general 
approach of crop acreage estimation and production. 

3.4.2 RESOLUTION FOR (S = 2, r = 2) ROTATION PATTERNS 

The rotation pattern (S = 2, r = 2) is presented In figure 3-1, which Is 
repeated as follows. 
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Figure 3-8.- 


Flow chart of a 


general estimation approach. 
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The design matrices of the related model, which is shown in equation (12), are 



or, with a = y /(1 + y) , 

“l a 0 0 0 o' 

a 1 0 0 0 0 

0 0 1 a 0 0 

H = (1 + y) 0 0 a 1 0 0 

0 0 0 0 1 a 

_0 0 0 0 a 1_ 




Then, 


(1 + Y)(l - 


-a 0 0 0 

10 0 0 
0 1 -a 0 

0 -a 1 0 

0 0 0 1 


TTTy 


1 + Y “Y 
-Y 1 + Y 


1 + Y 


1 + Y 0 

0 1 + Y 


1 + Y 


2(1 + y) “Y 

X*H”^X=y^^ -Y 2(1 + y) 


2(1 + y)J 


(X'H"lx)"^ = (1 + 2 y) 


2(1 + y) 


2(1 + y) “Y 
“Y 2(1 + y) 


The inverse matrix in the second member of the above equation is found by 
applying the result given in appendix B where 




a b 
b a 
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where, In this case, 
a = 2(1 + y) 


b a «Y 
n « 3 
and 


hence. 


that Is, 


= (a’^) 

^11 a (n ~ 2)b 

® ” La + (n ~ l)bj(a - b) 

.11 _ 2(1 + y) + (-y) 

® ^ it i 1 ■ "+' Y y + TF?rJi:^( 1 'TiTT yj 

2 + Y 

' wrrwi 


and 


,1j - -b 

‘ - ra'+ In - ni3-JTa-”¥r 


.1J 




So, 


(X'H“^X)"^ = 


1 + 2y 
2(r+ '3'yT 


"2 + Y 
Y 

- Y 


Y Y 
2 + Y Y 

Y 2 + Y. 


For this rotation pattern, the basic model 

a * Xa + Ub + le 


(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


can be written as 


L = a. + b + e. 
ts t s ts 


(45) 
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where 


t » 1, 2, 3 
s « 1, 2 , 3 
a * (aj^, ttg, ctg)' 

b - (bj^, bg, bg)’ 

That Is, 

^31 

Ai2 
^22 
A 23 

A 33 

U 

where 


1 0 0 
0 0 1 
1 0 0 
0 1 0 
0 1 0 
0 0 1 


«1 

“2 


1 0 0 
1 0 0 
0 1 0 
0 1 0 
0 0 1 
0 0 1 


bl 


eil 

®31 

®12 

®22 

®23 

®33 


(46) 


a = (A^, Agp ^^ 2 * ^ 22 * A 23 » A 33 ) 


(47) 


By step 4 of the acreage estimation In figure 3-8, the estimates of the A^g's 
are known as hence, the estimates of the a^'s are deduced from 
equation (28) which Is now 


(X'H"^X)'^X‘H"^a 


where 


and where 


0 / ist fV- f%J /v\ 

~ “ (All, Agp Aj^ 2 » A22. Agg, Aggj 


(48) 


(49) 


H = I + yUU', with Y = Og estimated from the pilot data, 
U Is given by equation (33), 

(X'H"^X)"^ Is given by equation (44), 
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X Is giveit by equation (32), and 
Is given by equation (34) with y. 


The variances of the estimates are deduced from equation (29) which is now 


Var(a) « (X'H“^X)“^Og = 

^ ® ® 2(2 + 3y) 


2 + Y 

A 

Y 

A 

L T 


A A p y A p 

where y ■ %/ %> estimated from the pilot data. 

This means 


'' (1 + 2 y )(2 + y ) 2 

var a. » s — 


2(2 + 3y) 


for t - 1, 2, 3, with T = 3. 


2 + Y 


2 + Y. 


(SO) 


(51) 


Emphatically, at the last year indexed t - T, 


Var(;^) . <M 2, i fH 2 ,^ ,,vl, ^ 


2(2 + 3y) 


With respect to this, note the following; 


(52) 


a. The datum T represents the total number of years utilized for the 

A 

multiyear estimator Oj. In the case of rotation with "return" 

(i.e., r < «), the maximum value of T is r + 1 since this corresponds to a 
complete cycle of the rotation. 

b. Remarkably, in the case of T = 3, if the rotation pattern cannot be as 
shown in figure 3-1, but instead will be shown in figure 3-9. 
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Then, the design matrices X and U will be as follows. 

’ 1 0 o” 

0 1 0 

0 0 1 
0 0 1 
_1 0 0 _ 

' 1 0 o" 
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which is the same as equation (36). Hence, equation (52) is obtained. 


-Y 

-y) 

2(1 + y) 


Var(L) = (L±M(2^) . 
' 2(2 + 3y) 


(56) 


(57) 


2(2 + 3y) 

However, the components of a°, which are given In equation (48), will be 
different In the case of figure 3-9, in which 


(58) 


0 rs* ^ \ I 

~ " \^ 11 * ^ 21 ’ ^ZZ* ^ 32 ’ ^ 33 ’ ^ 13 ) 


(59) 


rv fKt 




Components ^23’ ^31 longer needed, and and 

needed Instead, 


3.4.3 RESOLUTION FOR (S = 3, r = 2) ROTATION PATTERNS 

The rotation pattern (S = 3, r = 2) is presented in figure 3-2, which is 
repeated as follows. 
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(63) 

since, by the result given by appendix B, the inverses of the 3-by-3 matrices 
in equation (62) have their components 

Ji 


and 


(1 + y) + (3 - 2)v 

L (1 + y) + (3 - 1)yJL(1 + y) “ yJ 


_ 1 + 2y 
" 1 + 3y 


(64) 


So, 
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Vector a in the basic model, as shown in equation (12), for this rotation 
pattern is 


a = (A^j^, A 3 P A^p Aj^ 2 » ^ 22 ’ ^42’ ^13* ^23* ^33’ ^24’ ^34’ ^44^* 
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with the A^g's being the estimates of the A^s's of the above components, the 

a.'s of the estimates are deduced from 
t 

a = (X'H"^X)"^X'H"^aO (69) 


where 


a° = |a^, Ag^, A^^, Aj^ 2 » ^22’ ^42’ ^13* ^23’ ^33’ ^24’ ^34’ \4/' 

And, the variances of the estimates are deduced from the variance-covariance 
matrix 
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fs Aktp /AO , 

where y = <^ 5 /^ 6 ’ estimated from the pilot data. 
That means 


(71) 


VarlaJ ■ ‘1 " 3rl(V ... M 
3{3 + 8y) 


for t = 1, 2, 3, 4, and T = 4, or 
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(73) 


3.4.4 RESOLUTION FOR (S = 2, r = 3) ROTATION PATTERNS 

The rotation pattern (S = 2, r = 3) Is presented In figure 3-3, which is 
repeated as follows. 


3-23 















Then, with 


and 


Appendix C 
for T = 4, 


a = y/(1 + y), and (1 - a2)(l + y) * 


H = (1 + y) 


1 

a 


0 

0 

0 

0 


a 

1 


0 0 
0 0 


0 0 
0 
1 
a 


0 0 
0 0 


0 

0 0 
a 0 


0 

1 

a 


0 0 0 0 
0 0 0 0 


(1 + 2y)/(l + y), 

0 0 o' 

0 0 ' 0 

0 0 0 

0 0 0 

a 0 0 

1 0 0 

0 1 a 

0 a 1 



( 76 ) 


(77) 




1 + y 
1 + 2y 


2 -a 
-a 2 

0 -a 
-a 0 


0 -a 
-a 0 

2 -a 
-a 2 


“ i ' T ~ 2y C1rc(2,-a,0,-a) 


= Circ(2(l + Y),-Y»0,-y) 


(78) 


shows the elements of the inverse of the circulant in X'H"^X, which 
deduces 



( 79 ) 


3-25 



where y =■ ®b/^e’ estimated from the pilot data. 

3.4.5 RESOLUTION FOR (S = 3, r « 3) ROTATION PATTERNS 

The rotation pattern (S » 3, r » 2) is presented in figure 3-4, which is 
repeated as follows. 


Year number 



The design matrices related to the model shown in equation (12) are as follows. 




and 



and 


H*‘ ■■xrTyW 


I + a 
-a 
-a 
0 
0 
0 
0 
Q 
0 
0 
0 
0 
0 
0 
0 


-a 

I + a 
•a 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


•a 

I + a 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

I + a 
-a 
-a 
0 
0 
0 
0 
Q 
0 
0 
0 
0 


0 

0 

0 

I + a 
^a 
0 
0 
0 
0 
0 
0 
0 
0 
0 


Q 

0 

0 

-a 

-a 

I + a 
0 
0 
0 
0 
0 
0 
0 
Q 
0 


0 

0 

0 

0 

0 

0 

1 + a 
-a 
»a 
0 
0 
0 
0 
0 
0 


0 

0 

0 

0 

0 

0 

-a 

I ^ a 
-a 
0 
0 
0 
0 
0 
0 


0 

0 

0 

0 

0 

0 

-a 

-a 

I V a 
0 
0 
0 
0 
0 
0 


0 

0 

0 

a 

0 

0 

0 

0 

0 

l 4. a 
-a 
-a 
0 
0 
0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

-a 

I + a 
•a 
0 
0 
0 


0 

0 

0 

0 

Q 

0 

0 

0 

0 

•a 

^a 

I 4- a 

a 

0 

0 


0 

0 

0 

t 

0 

0 

0 

0 

0 

0 

Q 

0 

0 

I 4- a 

-a 

-a 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-a 

I + a 
-a 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

»ra 

-a 

I + a 


(83) 


Additionally, 


y'n"’^y ~ ^ 4* Y _ 
^ ^ “ 1 + 3 y 


1 


3(1 + a) 

-2a 

-a 


-a 

-2a 

-2a 

3(1 + a) 

-2a 


-a 

-a 

-a 

-2a 

3(1 + 

a) 

-2a 

-a 

-a 

-a 

-2a 


3(1 + a) 

-2-» 

-2a 

-a 

-a 


-2a 

3(1 + 

[1 + 2y) 

-2 

-1 


-1 


-2 

|(1 + 2y) 

-2 


-1 


-1 

-2 


2y) 

-2 


-1 

-1 

-2 


^(1 + 2y) 


-2 

-1 

-1 


-2 

^(1 

Y 


that is 

By appendix C, 


X'H"^X = 


1 + 3y 


C1rc(^l t 2it),-2,-l,-l,-2j 




If5 g-iY(il-l)2Tr/5 


-2 

-1 

-1 

-2 


(84) 


(85) 


( 86 ) 
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where 

Hence, wHh 


then 


k=l 


ad) =1(1 + 2 y) =1 + 


a(2) = -2 = a(5) 
a(3) = -1 = a{4) 


X- = X) a(k) = 3 /y 

^ k=l 


» |- + 6 - 2(.309 - ,9511) - (-.809 - .5881) - (-.809 + 5881) - 2(.309 + 9511) 

= I + 6 - 4(.309) + 2(.809) 

» ^ + 6.382 

* * 

Xj « a(l). + a(2)e"^^^^"^® + a(3)e"'''^’^^’'^® + a(4)e"^^’‘^’'^® + a(5)e”''®^^’’'^® 

»| + 6 - 2(-.809 - .5881) - (.309 + .9511) - (.309 - .9511) - 2(-.809 + .5881) 

a 1 + 6 + 4( .809) - 2( .309) 

= |-+ 8.618 


Similarly, 


X^ = - + 8.618 = X 3 


and 

Then, 


Xg = - + 6.382 = X 2 


, , , V,, -18x2it/5 -116 x2it/5 -124x2w/ 5 -132x2it/5 

* =) Hr * ’ ^-5T- * — "“T- 


1,1 . „( Y . -.809 t 5881 , .309 - .9511 . .309 + .9511 . -.809 - .5881 

i'v * 3" “3 ”3 1 — 3 "f 

^ i + g^glQ £+ g^g^g g ggg j 


i + 6.38Z 


(87) 


(88) 


(89) 


(90) 


(91) 


(92) 

(93) 
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That Is, 


“ i (7 * 3) 4 - 

t)5j & \Y / j J ^ g^3Q2 £ ^ g^g^Q j 


and 


(X'H-^x4-^1 + 3v) [1 H- rf ^TO y - 


1.618 I 


Then, for T = 5, 


Var a^, » ^(1 + 3 y) + 


0.618 


1.618 


3 + 8.618y 3+6 


IL^\ 

.382y/ 


A A p /a P 

where y = ^b/°e* estimated from the pilot data. 


(94) 


(95) 


(96) 


3.4.6 RESOLUTION FOR (S « 2, r = ~) ROTATION PATTERNS 
3. 4. 6.1 Case of T = 2 

The rotation pattern (S = 2, r = »; T’= 2) is presented below in figure 3-10. 

Year number 


Segment 

number 



1 

CM 

1 

X 


2 

X 

X 

3 


X 


Figure 3-10.- (S = 2, r = ») Patterns [2 retained 1; T = 2]. 
The design matrices related to the model shown in equation (12) are 

1 o' 

X = 1 1 0 

0 1 
0 1 


(97) 
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and 


U = 


1 

0 


0 0 
1 0 
0 I 0 

0 0 1 


(98) 


Then, 


H - I + yUU' - 


1 + Y 
0 
0 
0 


0 

1 + Y 
Y 
0 


0 

Y 

1 + Y 
0 


0 ■ 
0 
0 

1 + Y 


(99) 


and 


.-1 


1 


1 + 2y 


1 + 2y 

0 

0 

0 


(1 + y)'^ -y(1 + y) 

-y(i + y) (1 + y)^ 


0 

0 

0 

1 + 2y 


( 100 ) 


and 

X'H"^X = 
So, 


1 

(1 + y)(1 + W) 


(1 + 2y) + (1 + 2y)' 
“Y(1 ■*" y) 


-y(i + y) 

(1 + 2y) + (1 + y)^ 


( 101 ) 


(1 -I- y)(1 2y) 

C(l + 2y) + (1 + - y'^il - y)^ 


(1 + 2y) + (1 + Y.r 
y(1 + y) 


y(i + y) 

(1 + 2y) + (1 + y)^ 


( 102 ) 


Vector a in the basic model shown in equation (12) for this rotation pattern 
is 


£ = Aj^2> ^zz* ^23^' 


(103) 


With the A^g's being the estimates of the A^g's of the above components the 

(104) 


a^'s of the estimates are deduced from 


a = (X'H"^X)"^X'H"^aO 


3-31 


where, 


" (Ajj. Ai 2 . ? 22 , A 23 )' (105) 

Additionally, the variances of the estimates are deduced from the variance- 
covariance matrix 

Var(a) - (X'H"^X)"^ 0 g 


„ „ (1 Y)g + 2 y) ^ 

s* (J ii r»--iriM..-innni ^ ' r-- ..-rn r -■ ui ,1 

® C(l + 2 y) + U + yrr - Y^U + yr 


(i + 2 y) + (1 + y)‘ 
y(i + y) 


y(i + y) 

(1 + 2y) + (1 + y)^J 


A A A /a a 

where y « %/^q* estimated from the pilot data« 
That means 


Var 



(1 t y)(1 h- 2y)C(1 2y) + (1 + y)^l 

C(l + 2y) + (1 + yrr - y'^ii + yr 


for t s 1, 2, T = 2, or 


(106) 


Var 



(1 + y)(1 + 2 y)[(1 2 y) + (1 + y)^] 

^ ' ~ "A~~'h '"A aA ^ A 

[(1 + 2y) + (1 + y) ] - Y (1 + y) 


A A 

where y = 0 



as estimated from the pilot data. 


(107) 
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3. 4. 6, 2 Case of T « 3 

The rotation pattern (S ° 2, r a 3) is presented below In figure 3-11. 


Year number 



Figure 3-11.- (S « 2, r =. <») Patterns [2 retained 1; T = 3] 


The design matrices related to the model shown as equable '» (12) are 


and 


X a 


1 

1 

0 

0 

0 

0 


0 

0 

1 

1 

0 

0 


0 

0 

0 

0 

1 

1 



(108) 


(109) 
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I 


Then, 


and 


H » I + yUU' = 


i-1 = 


1 


1 + Y 
0 
0 
0 
0 
0 


1 + 2y 
1 + Y 
0 

0 

0 

0 

0 


0 

1 Y 

Y 

0 

0 

0 

0 

1 + Y 

-Y 

0 

0 

0 


0 

Y 

1 + Y 
0 
0 
0 

0 

-Y 

1 + Y 
0 
0 
0 


0 

0 

0 

1 + Y 
Y 
0 

0 

0 

0 

1 + Y 
“Y 
0 




F+Ty 


(1 + 2y) (1 + y)‘ 

1 + Y 

-Y 


0 


-Y 

2(1 + y) 

-Y 


So, 


(X'H'^X)“^ 


= (1 + 2y) 


(1 2y| 


+ (1 + y)' 


0 

0 

0 

Y 

1 + Y 
0 

0 

0 

0 

-Y 

1 + Y 
0 


0 

0 

0 

0 

0 

1 + Y 

0 

0 

0 

0 

0 

1 + 2y 
1 + Y 


( 110 ) 


( 111 ) 


-1 

(1 t 2t) t (1 t y)‘ 

1 + Y 


( 112 ) 


-1 


+ Y 
-Y 

0 


-Y 

2(1 + y) 


-Y 


(1 2y) + (1 + y)‘ 

1 + Y 


or, with H = (1 + 2 y) + (1 + y) > then 

. d/(l + y) 

(X'H"^X)'^ = (1 + 2y) 


-Y 

0 


-Y 

2(1 + y) 

-Y 


0 

-Y 

d/(l + y) 


(113) 
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which can be written as 


where 

and 


That means 


= (1 + 2y) 

lAl . 


f 1 = 1 + Y 

2 d(d - Y^) 


2 d - Y 

Y d 


Y d 
1 + Y 


JLl 


1 + Y (1 + y)2 ^ Y 


Y d 
1 + Y 


2 d - Y 


iu-lv\-l - (1 + y)(1 + 2y) 


(X'H"-"X) 


2 d(d - y)' 


2 d - Y 


r^Y 


1 + Y 
it + y) 


2 1 + Y 


rfvd 2d-Y' 


(114) 

(115) 


(116) 


(117) 


Vector a in the basic model shown In equation (12) for this rotation pattern 
Is 


a = (A^j^, A^o, Agg. A23, A23, A24)' 


(118) 


With the A^g's being the estimates of the A^g's of the above components (by 
step 4 of the flow chart In figure 3-10), the a^'s of the estimate are derived 
from 



a = (X'H’^X)‘’4'H"^aO 

(119) 

where 

“ (^11* ^12’ ^22’ ^23’ ^33’ ^ 34 )' 

(120) 


And the variances of the estimates are deduced from the variance-covariance 
matrix as follows. 
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Var a 




_ _2 (1 + y)(1 + 2y) 

- C 7 q A ' " " 

® Z did - Y ) 


2 d - Y 


1 + Y 


■'2 

Y 


where y 


"2 /"2 
" %/'^e’ 


1 + Y 
.2 




(1 + v)*' 1 + Y 


1 + Y 


-^d 2 d - y 


r.2 


as estimated from the pilot data, and 


:.x2 


d = (1 + 2 y) + (1 + y) 


That means, for T = 3, 


VarK] . al R ,?I 

® 2 d(d - Y^) 


A A 0 

where d = (1 + 2 y) + (1 + y)^ 

A AO 1^0 

and Y = %/ estimated from the pilot data. 

3.4.6 .3 Case of T - 4 

The rotation pattern (S = 2, r = ~; T = 4) is presented in figure 3- 
is repeated as follows. 


( 121 ) 


( 122 ) 


(123) 


5, which 
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Then, with a = y/U + y)» equation (125) can be written a-s 


hence. 


H = (1 + y) 


1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

a 

0 

0 

0 

0 

0 

0 

a 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

a 

0 

0 

0 

0 

0 

0 

a 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

a 

0 

0 

0 

0 

0 

0 

a 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 


■1 = 1 


1 + y 


1 

0 

0 


0 

1 


0 

-a 


1 - a^ 1 “ a^ 
-a 1 


1 - a^ 1 “ 


0 0 
0 0 


0 0 
0 0 


0 

0 

0 

0 

0 


0 

0 

0 

1 


0 

0 


-a 


l-a^ 
-a 1 


1 - a^ 1 “ a^ 


0 

0 

0 

0 

0 

1 


0 

0 

0 

0 

0 

-a 


( 126 ) 


0 

0 


0 

0 


1 - a^ 1 - a^ 
-a 1 


1 - a^ 1 - a^ 
0 0 


0 

0 

0 

0 

0 

0 

0 


(127) 
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and 


where 


X'H"^X 


1 

i + Y 


1 + - y — li-y 0 0 

1 - 1 - 



f 


1 

1 + Y 


g 

0 

0 


g 0 
2h g 
g • 2h 
0 g 


0 

0 

g 

f 


f = 1 + = C(i + 2 y) + (1 + y)^3/(i + 2y) 

1 - a^. 

g = -a/(l + a^) = - y (1 + y )/(1 + 2 y ) 

h = 1/(1 - a^) = (1 + y )^/(1 + 2 y ) 

By the Cramer rule in solving the system of equations, 

! fXj + gxg =0 

gx^ + 2hXg + gx2 = 0 

gxg + 2hX3 + gx^ = 0 

9 X 3 + fx^ - 1 


(128) 


the value will give the following value of the element (X'H"^X)"^44 


(X'H"4)"J = (1 + y) 


f 

9 

0 

1 

1 

0 

9 

2h 

9 

1 

1 

0 

0 

J— 

2h 

1 

_L 

0 

0 

0 

9 

! 

J_ 

1 

f 

9 

0 

1 

1 

1 

0 

9 

2h 

9 

1 

1 

0 

0 

9 

2h 

1 

1 

T' 

9 

0 

0 

9 

1 

1 

f 


which, by being applied to the calculation of the determinant of a partitioned 
matrix, is 


(X'H-^'XljJ = (1 


( 

+ Y)|f 


- (0 0 g) 


'f g O' 

-1 

fO]) 

9 2h g 



-0 g 2h, 


LgJ 


= (1 + y) ^ 
= (1 + y) 


f - 


2 fg^h - 

4 fh - fg '^ - 2 g^h 

4 fh - fg ^ - 2 g^h 

:2T“ 


r ' 

> 


+4f^^h - f^g^ - 4fg*"^h + g^ 


i .e., for T = 4, 


Varlouj.) = 0^(1 + ^)(4fh - fg^ - 2g^h)/(4f^h - f^g^ - 4fg^h + g^)| 


(129) 


where , 

Y “ 0 ^/ Og , as estimated from the pilot data 

f = C (1 + 2 ^) + (1 + y )^ 3/(1 + 2 ^) 
g = -y(1 + y)/(1 '+ 2y) 
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5 


S » (1 + y )^/(1 + 2 y ) 


3.4.7 RESOLUTION FOR (S = 3, r * “) ROTATION PATTERNS 
3. 4. 7.1 Case of T = 3 

The rotation pattern (S = 3, r = T = 3) Is presented below in figure 3-12. 

Year number 


Segment 

number 


X 

1 

2 

CO 

1 

X 



2 

X 

X 


3 

X 

X 

X 

4 


X 

X 

5 



X 


Figure 3-12 - (S « 3; r * «) Patterns [3 retained 2; T = 3]. 
While this pattern shows the notation for only three years, 
the related pattern given in figure 3-6 shows it for four. 


The design matrices related to the basic model, shown in equation 12, are 

‘l 
1 
0 


1 

0 

0 

0 

0 

0 


0 

0 

1 

0 

1 

0 

1 

0 

0 


0 

0 

0 

0 

0 

1 

0 

1 

1 


(130) 
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and 


U = 


1 

0 

0 

0 

0 

0 

0 

0 

0 


0 0 0 0 
10 0 0 
10 0 0 
0 10 0 
0 10 0 
0 10 0 
0 0 10 
0 0 10 
0 0 0 1 


Then, with a 


= y(1 + y)» u 


H = (1 + y) 


= (1 

~1 

0 

0 

0 

0 

0 

0 

0 

0 


- a)(l + 2a) = (1 + 3y)/(1 + y)^ 

0000000 o” 

1 a 0 0 0 0 0 0 

a 1 0 0 0 0 0 0 

OOlaaOOO 
0 0 a 1 a 0 0 0 

OOaalOO O 
0 0 0 0 0 1 a 0 

0 0 0 0 0 a 1 0 

0 0 0 0 0 0 0 1 


(131) 


(132) 



0 

0 



H-1 - 1 

H -r+Y 


0 • 0 
0 0 

0 0 

0 0 


0 


-a 



0 

0 

0 

0 

0 

0 


0 0 

0 0 

0 0 

1 + a -a 
u u 

zl 1 + a 
u u 

zi zi 

u u 

0 0 

0 0 

0 0 


0 0 0 0 
0 0 0 0 

0 0 0 0 


-a 


u 


0 


0 


-a 

u 

1 + a 
u 

0 

0 




0 0 0 


0 

0 

0 

0 

0 

0 


(133) 
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So, 


X'H"^X = — ^ 


r+Y 


1 + 1 « 


1 - a 

j, a 

u " 




1 - a‘ 

-a 
“iT 


which can be written 


^ 1 


-a 

u 


1 + a ~a a 




1 - a 
-a a 


1 - a 


u 


1 - a 


„ 1 a. 1 X 1 + a 

7 1 + - -J + -— 


1 - a 


(134) 


X I u'"^V - 1 

I* ^ " 1 + Y 


"A 

B 

D" 

B 

C 

B 

.0 

B 

A. 


(135) 


where 


1 *“ 3 


n _ -a a _ -y(1 + y)[(1 + 2y) + (1 + 3y)] 

(T+aY)(l + 3 y) 


(137) 


P - o 1 . 1 + a _ (1 + y)C2(1 + y)(1 + 3y) + (1 + 2y)^1 

^ = u (1 + 2 y)(1 + 3y) 

i “ a 


and 


n - -a _ ^-<(1 ’H y)(1 + 2y) 
^ “ u " (1 + 2 y)(1 + 3y) 


(139) 


In applying the definition of the Inversion of a matrix, the three diagonal 
elements of the Inverse of X'H“^X will be 


AC - B 


(X'H"^X):} = (X'H"^X):J = (1 + y) o 

(A - D)C(A + D)C - 2B^] 


(140) 


and 


(X'H"^X)22 = (1 + y) 


(141) 


A + D - B‘ 
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That means 


Var * (1 + y)*ii 

• / A 


/OO 0 ?, 
(AC - B^) 
^'OTTT 


07770 07D 

(A - D)C(A + D)C - 2B‘^] 


for T * 3 , where 

A A p /h p 

Y * %/ 3s estimated from pilot data 

A » (1 ^ 2y)(1 •<• 3y) (1 + y)^(1 3y) + (1 + y)(1 + 2y)^ 

(1 + 2y)(1 + 3y) 

B = -y( 1 + v)C(l 2 y) + (1 3y)3 

(1 + 2y)(1 + 3y) 

C = (1 + Y)t2(i + y)(1 + 3y) + (!•<• 2y)^] 

(1 + 2^)(1 + 3 ^) 

and 


(142) 


(143) 

(144) 

(145) 


D = jHl . t 2 y) . (14g) 

(1 + 2y)(1 + 3y) 

Note: In applying the Cramer approach to directly solve for X3 In the system of 
three equations 

! AXj^ + BXg + DXg = 0 

BXj^ + CXg + BX3 = 0 (147) 

DXj^ + BXg + AX3 = 1 

the value of the element (X'H’"^X)23 will be as follows. 
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(X'H"^X)33 * U + y) 


A 

B I 
1 

0 

B 

mmrn 

iJ 

0 

D 

B ! 

1 

A 

TT 

D 

B 

c 1 

B 

D 

B { 

A 


and by calculating the determinant of a matrix that Is partitioned, 
that Is 


(X'H’^X)'^ - (1 + y) 



DMI-1 
B 


A - (0B)(* ®)-l( 

2 

(X'H'^Xlgi « (1 + y) y - 5 ^ 

1A(AC - B ) - DT + 2B D - AB'^I 


- n + -1 - B 

- U + YJ "7J 5 7 

ii u/> r> A n M T\^i 


AC B 
2 


AT - 2AB‘- - DT + 2BT 
2 


— /I j. \ AC “ B 

= 11 + y) — 5 -T 7 

/ A ^ on ^ 


= (1 + y) 


(A - D^)C - 2B^(A - D) 
AC - B^ 


(A - D)C(A + D)C - 2B^] 
which Is the same as equation (140), 


(148) 


(149) 


(150) 


3. 4. 7. 2 Case of T - 4 

The rotation pattern (S = 3, r = «; T = 4) Is presented in figure 3-6, which 1s 
repeated as follows. 


Segment 

number 


Year number 


X 

1 

2 

3 
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X 


(S a 3, r ** «) Patterns C3 retained 2} T => 4] 


The design matrices related to the basic model, as shown on equation (12), are 


X = 


1 

1 

0 

1 

0 

0 
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0 

0 

0 


0 

0 
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0 

1 

0 

1 

0 

0 
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0 

0 


0 
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0 

0 

1 

0 

1 

0 

1 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

1 


(151) 
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10 0 0 
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0 0 0 0 
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0 0 1 
0 0 1 
0 0 1 
0 0 0 


0 0 
1 0 


1 0 


0 0 10 0 


0 0 0 0 1 0 

0 0 0 0 1 0 

0 0 0 0 0 1 ^ 

Then, with * y/(1 ■** y)> u » (1 - a)(l + 2a) ■ (1 + 3y)/(1 + y)^! 

100000000000 
OlaOOOOOOOOO 
OalO '0 0000000 
OOOlaaOOOOOO 
OOOalaOOOOOO 
1 -|»^+y) OOOaalOOOOOO 
OOOOOOlaaOOO 
OOOOOOalaOOO 
OOOOOOaalOOO 
OOOOOOOOOlaO 
000000000 a 10 
000000000001 





where A, B, 
E 


0 


are 

1 

1 - 


defined in equations (137), (138), and (139), and 


3-48 


Therefore, 


(X'H"^X)4J 


* (1 + y) 


Hence, for T = 4, 


A 

B 

D 1 0 

B 

E 

2D 1 0 
1 

D 

2D 

j_L“_ 

0 

0 

B ll 

A 

B 

D 1 0 

B 

E 

2D 1 D 

D 

2D 

E 1 B 

0 

D 

B 1 A 


(1 + y) <A - 


( 

“A B D1 

-1 

’Ol) 

{A - (0 D B) 

B E 2D 


D } 

( 

D 2D E 


Bj) 


n\-l 


(156) 


A A. { 0 0,0*10 1 _1 2 

Var = (1 + y) |A - b'M ^b| 


(157) 


where 


0 

A 


m = (1 + 3 y)/(1 + y) 

[(1 -I- 2y)(1 + 3y) (1 + y)^(1 + 3y) + (1 + y)(1 + 2y)^] 

(1 + 2^)(l + 3^) 


(158) 


B = -y(1 + y)[( 1 + 2y) + (1 + 3t)] 

(1 + 2^)(1 + 3 ^) 


(159) 


D - -y(1 y)(1 + 2y) 

(1 + 2^)(1 + 3y) 


(160) 


E = (1 + y)C(1 + y)(1 + 3y) + 2(1 + 2y) ] 
(1 + 2^)(1 + 3y) 


(161) 
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(162) 


and 




0 0 
D 2D 



0 0 0 0 
b' = CO D B] 

r\j 


(163) 


3.4.8 RESOLUTION FOR (S = 4, r = ~) ROTATION PATTERNS 


3. 4. 8.1 Case of T = 3 


The rotation pattern (S = 4, r = T = 3) is presented below in figure 3-13. 


Year number 


Segment 

number 



Figure 3-13.- (3 = 4, r = «) Pattern [4 retained 3; T = 3]. 
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Then, with a = y/(1 + 

y)» 
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H = (1 + y) 
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0 
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1 _ 

and by applying the results given ’ 

in 

appendix 

A with 







U 

= (1 

- 

a)(l 

+ 2a) 







= (1 + 3y)/(1 + y)^ 

then 


(166) 



(167) 
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and 


1 + 2 — jj— 




1 - a‘ 


1 

T-fJ 


1 - a 
-2 

F G 2D 
G H G 
2D G F 


7“ ‘ u 


2 4 2(- 


a o a 

IT7 “ 

I X 1 + 3a 


1 - a 


7 u 


-2i. 

u 


a -2a 


1 - a 


7 •• u 


a . 2 4 1 + —i— y + 2 Ltl 


1 - a 


7’ ‘ u 


1 - a‘ 


(168) 


where 


(1 + 2y)(1 + 3y) + I 

[1 + y)^i 

fl + 3 y] 

» + 2(1 + y)(1 + 2y)2 

< 

;i + 2 y)i 

[1 + 3y1 

1 


G = - 


H 5 2 


a o a _ -y(1 + 2y)[2(1 + 2y) + (1 + 3y)] 

- 2 ir -Trrr+ -2iT (i—IyITi'+'^y) 

(r^ 


1 + a 


, X . - \ _ 2(1 + y)C(1 + y)(1 3y) + (1 + Z-^ri 

J + ■(r-.-^JTr TTalj n-x x ''^.v »— 


(1 + 2y)(1 + 3y) 


D = - J 1 as in equation (140), 

Therefore, 


(171) 

(172) 


=« (1 + y) 


F 

g! 

o| 

G 

HI 0 


—J 

2D 

6! 1 

F 

G j 2D 

6 

Hj G 

2D 

g| 

F 


(1 + y) jiF - (2D G) §) 


FGW (2D), -11 


.. . , (FH - G'^) 

a (1 + Y) ■ ' ? ' ? 0 2 

F(FH - 6^) - 4D*^H + 40G^ - FG^ 


(1 + y) 


FH - G‘ 


« (F - 2D)[(F + 2D)H - 2G^] 


(173) 
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In fact, 




Hence, for T » 3, 



( 174 ) 


where 



Y = %/^e* estimated from the pilot data 


» 

\ 

^ _ [(1 + 2y)(1 + 3y) + (1 + y)^(1 + 3y) + 2(1 + y)(1 + 2y)^] 

(175) 

\ 

1 

(1 + 2^)(1 + 3^) 


J _ -y(1 + 2^)[2(1 + Zy) + (1 + 3^)] 
(1 + 2^)(1 + 3^) 

(176) 


^ _ 2(1 + y)C(1 + y)(1 + 3y) + (1 + 2y)^] 
(1 + 2?)(1 + 3y) 

(177) 


g _ -y(1 + y)( 1 + 2^) 
(1 + 2^)(l + 3^) 

(178) 


3.4. 8. 2 Case of T = 4 

The rotation pattern (S = 4, r = «; T = 4) is presented •'’n figure 3-7, which 
is repeated below. 
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Figure 3-7.- (S * 4, r « ») Patterns [4 retained 3]. 
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Then, with a = y/1 + y» 


H - (1 +y) 


1 

0 
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0 

0 

0 

0 
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0 

0 
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0 0 1 a a 0 
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0 0 a a 1, 0 

0 0 0 0 0 1 

0 0 0 0 0 a 
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0 , 0 0 0 0 0 
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0 0 0 0 0 
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0 0 0 0 0 
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a a 0 0 0 
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0 0 a 1 0 
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(181) 


and, by applying the results given In appendix A, we denote equations (182) 
through (184) 

n = (1 - a)(l + 3a) = (1 + 4 y)/( 1 + 
u = (1 - a)(l + 2a) = (1 + 3 y! '1 + y)^ 

where 



J = 



1 - a . 1 2a _ 
u n “ 


1 + 


(1 + 
TT 



1 + Y . (1 + y) (1 + 3y) 
1 + 3'y 1 + 4y 


(185) 


a a a ^ yU + y) y(1 + y) y(1 + y) 

~ 2" “ u “ n " U + <^yT " 1 + 3y 1 + 4y 

1 •• d 


(186) 


a a _ y( 1 t_ll y(1 y) 

u ” n ~ ” 1 + ”3y 1 + 4y 


(187)' 


a _ y(1 + y) 

- - ^ - 1 + 4y 


(188) 


«a 
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( 189 ) 


n «a a _ 2 y(1 + y) y(1 + y) 

Q « -2-jj - - = - nrTlY”* - TTTT 


. _ 1 . „ 1 + a . 1 + 2a _ (1 + y) . 2(1 + y)(1 + 2y) , (1 + 

p « ^ + 2 — » -- rHr ^ — r+-TY * ^ 

1 - a 


Therefore, 


(X'H“^X)^J = (1 + y) 


J K L I 0 

K P Q I 0 
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R L K I 1 
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l; R 


J K 

K P Q I L 

L Q P I K 

R L K ! J 
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■J 
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-1 

r«i)‘ 

= (1 t y) < 

'o - (R L K) 

K 

P 

Q 



i 

1 

.C 

Q 

P- 


UJ) 


Hence, for T = 4, 


A A 

Var Oj = (1 + y) 



•1 ol“l 
VI 


where 


Y = %f estimated from the pilot data 

°i - 1 J. (1 y)^ j. 1 y + (1 + y)(1 + 3y) 


y)(1 + 3y) 
+ 4y 

(190) 


(191) 


(192) 





( 194 ) 


Kb- + y) _ y( 1 ± A _ y(1 + y) 

1 + 2y 1 + 3y 1 + 4y 

R a “Yd + y)/(1 + 4y) 

L a “y(1 + y)/(1 + 3y) - y(1 y)/(1 + 4y) 

Q a “2 y(1 + y)/(1 + 3y) - y(1 + y)/(1 + 4y) 


(195) 

(196) 

(197) 


P a (1 + ^)^/(l + 2 y) + 2(1 + y)(1 + 2^)/d + 3?) + (1 + y)(1 + 3y)/(1 + 4y) 


and 


0 , ,0 0 0 , 

V* a R L K) 

rw 


0 0 0 
J K L 

0 0 J 

K P Q 


(198) 

(199) 


.0 0 0 , 

LL Q PJ 


( 200 ) 


3.4.9 RESOLUTIONS FOR OTHER SPECIFIED ROTATION PATTERNS 


Following are some other specified rotation patterns of S = 4 segments per 
year, the results of which are to be presented and applied In the cases of 
necessary collapsions mentioned In section 3.3.2. 


3. 4. 9.1 Resolution for (S =4, r = 3) Rotation Patterns 

The rotation pattern (S = 4, r = 3; [4 retained 3]) Is presented In 
figure 3-14. 
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Figure 3-14.- (S = 4, r = 3) Patterns [4 retained 3]. 


Some results are 


X'H 


Iv > 1 


1 + 4y 


C1rcC4(l + 3y), -3y, -2y, -2y, -2y, -3y] 


(X 


„-lvrl - 1 + 4y /I . 1 

» X'66 ^ 4 + i6y " 


/J 


4 + 14y - /3y 4 


— ) 

+ lOy + 5 /3y/ 


and, for T = 6, 


= 1 + 4y /I ^ __1 /J /J \ 2 

r 4 + 16y 4 + 16y - /3y 4 + 10^ + 5 /3y / ® 

where y = %/°e’ estimated from the pilot data. 

3. 4. 9. 2 Resolution for (S = 4, r = 4) Rotation Patterns 

The rotation pattern (S = 4, r = 4; [4 retained 3]) is presented in 

figure 3-15. 
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Figure 3-15.- (S = 4, r = 4) Patterns [4 retained 3]. 


Some results are 


«T-i-i:-Circ[4(i + 3y). - 3y, -2 y, -y, -y, -2y,-2y, -3y3 (204) 


/V I u“l •• 1 _ 1 + 4 y /I X 1.470 

(X H X)^7 7 — 4 +15.3' 


1.802 


0.445 


(205) 


and, for T = 7, 


" .. 1 + 4y/1 X 1.470 1.802 0.445 

Var «-|- 7 I -r + - 7T • X - 

' \ 4 + 15.357Y 4 + 15.602Y 4 + 10.951y 


(206) 


'2 rz 


where y = ^b/'^e’ estimated from the pilot data. 


3. 4. 9. 3 Resolution for (S = 4, r = 2) Rotation Patterns 


The rotation pattern (S = 4, r = 2; [4 retained 2]) is presented in 
figure 3-16. 





























Year number 



Figure 3-16.- (S = 4, r = 2) Patterns [4 


retained 2]. 


Some results are 


X'H"^X 


2 

TTTy 


.-1 


2(1 + y) -Y -Y 

-Y 2(1 + y) -Y 
-Y -Y 2(1 + y) 


= 2(X'H-‘X)p,g_ 3.1 

= l(X'H'^X)p5g 3_i 


and, for T = 3, 


Var 


oc^i =7Varia^ip.g^ 1 


_ (1 + 2 y )(2 + y) 2 


4(2 + 3y) 


(207) 

(208) 


(209) 


3. 4. 9. 4 Resolution for (S = 4, r = 3) [4 retained 2] Rotation Patterns 

The rotation pattern (S = 4, r = 3; [4 retained 2]) is presented in 
figure 3-17. 
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Figure 3-17,- (S = 4, r = 3) Pattern [4 retained 2], 


Some results are 

X'H"^X = C1rcC2(l + y), -Y, 0, -y] 


(X'H“^X)"^ = 3 ( 


and, for T = 4, 



3. 4. 9.5 Resolution for (S = 4, r = ”) Rotation Patterns [4 Retained 2] 

The rotation pattern (S = 4, r = «; T = 3, [4 retained 2]) is presented in 
figure 3-18. 


3 - 6 ' 



Year number 
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number 



Figure 3-18.- (S = 4, r = «) Patterns [4 retained 2; T = 3]. 


Some results are . 

f H 1 + 1/(1 - a^) = C(1 + 2y) + (1 + y)^]/(1 + 2y) 
g = -a/(l - a^) = -y(1 + y)/(1 + 2y) 
h = 1/(1 - a^) = (1 + y)^/(1 + 2y) , as in equation (128) 

y g 0 


x'H-‘x =rlr 


(X'H"^X)’^ = -^-5^ 


2f(fh - gn 


9 2h g 
Lo g fj 


2fh - -gf 

.2 


-gf 

2 


-gf 


-gf 2fh - g^ 


and, for T = 3, 


. » ? , -'. , 00 . 0 ? 0 00 op 

Varlop) - cTg (1 + Y)(2fh - g"^)/4f(fh - g'^) 


(213) 


(214) 


( 215 ) 
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where 


f = [(1 + Zy) + (1 + ^)^]/(l + 2^)\ 

I = -;(1 + ;)/(! + 2 ;) ( ( 216 ) 

h = (1 + ^)2/(l + 2^) ) 

3.5 EMPIRICAL OPTIMIZATION,' WITH PRECISION ACHIEVABLE BY ROTATION SAMPLE 

design EstiWATIDN 

As frequently mentioned In previous sections, the estimated values y of y are 
obtained from Lapdsat (pilot) data. In all numerical evaluations of the 

A 

variance formulas, these estimated values y are to be substituted. 

In addition, as discussed in section 1.5.1, the objective of rotation sample 
designs is to utilize the "consistency" of the wheat acreage of a particular 
segment from year to year in order to improve the accuracy of the current 
year's wheat acreage estimate. This accuracy may be seen from the gain in 
precision which is attainable compared with the aggregation based only on the 
current year's estimates. 

The statistic presenting the gain in precision is named the variance reduction 
ratio R. 

3.5.1 VARIANCE-REDUCTION RATIO R 

3. 5. 1.1 Variance v(Tj) of Estimators Based Only on Current Year Acquisitions 

For any value of T, for example T =3, the estimator can be shown as 
follows. 




■3-66 



Yj = •j(0...0|0...0|l...l|0...0|0...0) 


A 


11 


IS 


T1 


TS 


T h'a 
s ~ ~ 


= S A.„/S 


J=1 


So, 


Var^Tj.) = 311 ' Var(a)h 


By equation (13), Var(a) = la^ + tf^UU' , equation (218) can be written as 


Var(Tp) = ^ h'h Og + h'UU'ba| 


S $“ 

where U is any design matrix related to equation (218) and 

h'h = S 


h'UU'h = S 


That means 


Var(T,.) = 4/5 + 4 /s 

“ ^ 4/4) A 


(217) 

(218) 

(219) 

( 220 ) 


( 221 ) 
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( 222 ) 


2/2 

and with y h %/^q* 


Var|Tjj = a^d + y)/S 


3. 5. 1.2 Definition of the Variance-Reduction Ratio R 

In order to compare the variance Var|ayj of the rotation sample current-year 
estimator oy to the variance Var|yjj| of the only-based current-year estimator 
Tj, the variance-reduction ratio R is Introduced as follows, 

R 5 Var|otyj/Var|Tyj 

That is, R = Var(iy] | (1 + Y)<Je/s| (223) 

This ratio R, then, represents the gain of precision achievable by the 
rotation sample design estimation; and, Is the variance reduction ratio for 
the rotation pattern according to any figure 1, 1 = 1, 2, •»*, 12. 

3. 5. 1.3 Analytic Formulae for R 

A 

Based on the formulae of Var Oj given in section 3.5, the formulae shown as 
equation (223) provides the following various analytic formulae for R 
according to different specified rotation patterns which are noted for 
figures 5-9 through 5-18 as follows. 

a. Figure 3-1. (S = 2, r = 2; T = 3, [2 retained 1]) 

aid + 2^)(2 + ^)/2(2 + 3d 

R =_e 

^ Og(l + y )/2 


by equation (52). 
That Is, 


Rj = (1 + 2d(2 + d/(l + d(2 + 3d 


(224) 
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Figure 3-2. (S = 3, r = 2; T = 4, [3 retained In 2]) 


clil + 3 y )(3 + 2 y )/ 3{3 + 8 y ) 

D =-S , 

Ogd + y)/3 

by equation (74). 

That is, 

Rg “ (1 + 3 y )(3 + 2 y )/(1 + y )(3 + 8 y ) 

c* Figure 3-3. (S = 2, r = 3; T = 4, [2 retained 1]). 

ahl + 1/(1 + 2 y ) - 2/(1 + y )](1 +■ 2 y )/8 
^ Ogd + y )/2 

by equation (79). 

That is, 

R 3 = Cl + 1/(1 + 2 ^) - 2/(1 + ^)](1 + 2?)/4(l + ^) 


(225) 


(226) 


Figure 3-4. (S = 3, r = 3; T = 5, [3 retained 2]). 

o^Cl/3 + .618/(3 + 8.618^) - 1.618/(3 + 6.382y)](1 + 3^)/5 
R^ “ o ^ 

Ogd + y)/3 

by equation (96). 

That is. 


R 4 = 3C1/3 + .618/(3 + 8.618^) - 1.618/(3 + 6.382y)](1 + 3^)/5(l + y) 


(227) 
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e. Figure 3-10. (S = 2, r = ~; T = 2, [2 retained 1]). 

By equation (107), 

„ _ 2(1 + 2 ^) C (1 + 2 ^) + (1 + 

C(1 + 2y) + (1 + yV-1 - Y^d + y)^ 

f. Figure 3-11. (S = 2, r = «; T = 3, [2 retained 1]). 

By equation (123), 


R = (1 f 2 y) {2C(1 + 2y) + (1 y)^] - Y^> 

C(1 + ?.y) + (1 + y)^]{[(1 + 2^) + (1 + ^)^] - 

9* Figure 3-5 (S = 2, r = ~; T = 4, [2 retained 1]) 

By equation (129), 



00 002 020. 
2(4fh - fg^ - 2g‘^h) 

" ()p o (jwir ojr 

4f^h - - 4fg^h + g^ 


where 

?«■[(! + 2^) + (1 + yfy/il + 2y) 
g = - y (1 + y )/(1 + 2 y ) 
h = (1 + y )^/(1 •’* 2 y ) 


h. Figure 3-12. (S = 3, r = <»; T = 3 [3 retained 2]). 
By equation (142), 


R 


6a 


00 02. 

3(AC - B'^) 

0 0 0 0 "0 07 

(A - D)C(A + D)C - 28*^ 


*i. Figure 3-6. (S = 3, r = T = 4, [3 retained 2]). 


( 228 ) 


(229) 


(230) 


( 231 ) 


By equation (157), 


. 1 . Figure 3-13. (S = 4, r = »; T = 3, [4 retained 3]) 


By equation (174), 


R = 4-CFI1 - G^ ) 

(F - 2D)[(F + 2D)H - 2G'^] 


( 233 ) 


0 0 0 0 

where F, G, H, and D are defined in equation (175) through (178). 
k. Figure 3-7, (S = 4, r = «; T = 4, [4 retained 3]. 

By equation (192), 


1 0 0, 0_1 0 ) _i 

R 7 = 4 jj - v' N ^ VI 


(234) 


0 0 

where J, v, N are defined In equations (194) through (201). 
1 . Figure 3-14 . (S = 4, r = 3; T = 6 , [4 retained 3]). 

By equation (203), 


Z3- 


/r 


R - 2(1+4y) /1,. 1 ^ 

^ 3(1 + y) 4 + 16y 4 + 14y ~ Y 4 + IOy + 5 /T y> 


(235) 


’’’• Figure 3-15. (S = 4, r = 4; T = 7, [4 retained 3]). 
By equation (206), 

1.802 


0.445 


R = 4(1 + 4y) n + 1.470 _ 

^ 7(1 + y) 4 + 15.357Y 4 + 15.602y 4 + 10.951y 

Figure 3-16. (S = 4, r = 2; T = 3, [4 retained 2]). 

By equation (209), 

Rj^q = (1 + 2 y )(2 + y )/(1 + y )(2 + 3y) = Rj^ 

0 . Figure 3-17. (S = 4, r = 3; T = 4, [4 retained 2]). 

By equation (212), 


(236) 


(237) 


Rl, = /l ^ -i-. - 

4(1 + y) \ 1 + 2y 1 + Y, 


R. 


(238) 
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P* Figure 3-18. (S * 4, r * «; T = 3, [4 retained 2]). 

By equation (215), 

,00 Oo 0 00 Qo 

R ^2 = (2fh - g^)/f(fh - g^) (239) 

0 0 0 

where f, h, and g are defined in equation (230), 

3.5.2 VALUES OF VAR|o^j AND R FOR A SPECIFIED VALUE OF y 

Table 3-1 provides a summary of the rotation pattern figures given in 
section 3 of this document. In addition, tables 3-2, 3-3, 3-4, and 3-5 show 
the values of Var^Ojj and R corresponding to the first seven rotation designs 
for two specified values of y, Y = 2.01 and y = 3.38. 

Based on a limited data bank, the estimation of y was performed in two cases 
as follows. 

a. Case 1 

A 4-year set .of data in which all acquisitions for all years were used; for 

# A 

’this case, the estimate was y = 2.01. 

b. Case 2 

A situation such as in Case 1, but one in which the first acquisition estimates 
in the fourth year were used; for this case, the estimate was y = 3.38. 

At this point, it is necessary to note that the estimate of y computed is not 

^2 ^2 

strictly a ratio of crj^ and which is computed for a single stratum. The 

data are not adequate for such an analysis. Instead, the available data for 

2 2 

all strata were pooled and a simulation of and was computed by fitting 
the model 

^tsc “ “t * '’s ®tsc '240) 

where the subscript c denotes the stratum (CRD) of the segment, the a^, y., 

u C 

and (ay)^(. are year, stratum (CRD), and interaction constants, respectively, 
and the bg, and e.j. 5 ^ are, as before, segment and error variables. 
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TABLE 3-1.- SUMMARY OF ROTATION' PAHERN FIGURES 
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TABLE 3-2.- VALUES OF VAR(i^) FOR ^ = 2.01 


No. of years 
1n rotation, 
T 

Values for figures 
with return 

Values for figures 
without return 

3-1 

3-2 

3-3 

3-4 

3-5 

3-6 

3-7 

2 

1.32 

0.90 

1.32 

0.90 

1.32 

0.90 

0.69 

3 

1.25 

.87 

1.29 

.83 

1.29 

.83 

.64 

4 


.86 

1.17 

.80 

1.285 

.82 

.61 

5 




.78 

1.284 

.81 

.59 


TABLE 3-3.- VALUES OF VAR(aj.)FOR y = 3.38 


No. of years 
In rotation, 
T 

Values for figures 
with return 

Values for figures 
without return 

3-1 

3-2 

3-3 

3-4 

3-5 

3-6 

3-7 

2 

1.81 

1.25 

1.81 

1.25 

1.81. 

1.25 

0.97 

3 

1.72 

1.22 

1.73 

1.13 

1.73 

1.13 

.88 

4 


1.21 

1.54 

1.0826 

1.71 

1.0833 

.817 

5 




1.056 

1.70 

1.060 

.787 
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TABLE 3-4 VALUES OF R FOR y “ 2.01 


No. of years 
1n rotation 
T 

Values for figures 
with return 

Values for figures 
without return 

3-1 

3-2 

3-3 

3-4 

3-5 

3-6 

3-7 

2 

0.877 

0.897 

0.877 


0.877 

0.897 

0.917 

3 

.831 

.867 

.857 

.827 

.856 

.827 

.850 

4 


.857 

.777 

.797 

.854 

.817 

.811 

5 




.777 

.853 

.807 

.784 


TABLE 3-5.- VALUES OF R FOR y = 3,38 


No. of years 
in rot at i on i 
T 

Values for figures 
with return 

Values for figures 
without return 

3-1 

3-2 

3-3 

3-4 

3-5 

3-6 

3-7 

2 

0.826 

0.856 

0.826 

0.856 

0.826 

0.856 

0.886 

3 

.785 

.836 

.790 

.774 

.790 

.774 

.804 

4 


.829 

.703 

.741 

.781 

.742 

.746 

5 




.723 

.776 

.726 

.719 



































In reviewing table 3-2 through 3-5, the following note iis pertinent. As 
mentioned in section (3, 5. 1.2. 2), in the case of rotation with return 
(i.e., r < “>), the maximum value of T is r + 1 since this corresponds to a 
complete cycle of the rotation; T is the total number of years utilized for 
the rotation sample (multiyear) estimator. No entry Is shown for T = 1 since 
this would be Identical with the single-year estimator. 


3.5.3 COMPARISONS BETWEEN ROTATION PATTERNS 

Based on the values of R, instead of the values of Var Oj., either table 3-4 or 
table 3-5 will first reveal the difference between rotation patterns with 
return and without return. In order to have a fair comparison, the values of 
R in th 6 se two tables will be presented in tables 3-6 and 3-7, the lines of 
which show the patterns with the same number of observations. 


Once the rotation patterns with return (i.e. r < ») reach their cycles, their 
variance reduction ratios have values which are less than those of the 
rotation patterns without return (i.e., r = »). 


In detail, tabl 
R]^ = ,'i31 < R 5 

R 3 = .777 < Rg 

h = < "6 

^2 ~ ^ 


es 3-6 and 3-7 indicate the following in the case of y = 2.01. 
= .856 for T = 3, S = 2, 6 observations in total 

= .854 for T = 4, S “ 2, 8 observations in total 

= .807 for T = 5, S = 3, 15 observations in total 

= .817 for T = 4, S = 3, 12 observations in total 


In the case of y 
Rj^ = .785 < R 5 = 
R^ ~ .703 ^ Rg — 
R 4 = .723 < Rg = 
R 2 ~ .829 ^ Rg — 


= 3.38, tables 3-6 and 3-7 indicate the following. 
.790, for T = 3, S = 2, 6 observations in total 

.781, for T = 4, S = 2, 8 observations in total 

.726, for T = 5, S = 3, 15 observations in total 

.742, for T = 4, S = 3, 12 observations in total 
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TABLE 3-6.- VALUES OF R IN ROTATION PATTERNS 
WITH AND WITHOUT RETURN FOR y = 2.01 


1 



0.877 

.856 

0.877 

.856 

.854 

0.897 

.827 

.817 
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TABLE 3-7.- VALUES OF R IN ROTATION PATTERNS 
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Hence, when they reach full cycle, the rotation patterns with return are 
better than those without return. Therefore, an optimal rotation design Is 
one which in most situations has the rotation pattern with reti rn and with 
minimum value of R. 

3.5.4 THE OPTIMAL ROTATION PATTERN 

The optimal rotation pattern will be a rotation pattern with return and with 
minimum value of R over all situations. 

A 

3. 5.4.1 Values of R for Rotation Patterns with Return for y » 2.01 and 3.38 

In the case of specified values y = 2.01 and y = 3.38, the values of R among 
rotation patterns with return are presented in table 3-8. 

TABLE 3-8.- VALUES OF R IN ROTATION PATTERNS 
WITH RETURN FOR Y = 2.01 and 3.38 


. No. of years 
in rotation, 
T 

R values per figure, 
where y = 2.01 

R values per figure, 
where y = 3*38 

3-1 

3-2 

3-3 

3-4 

3-1 

3-2 

3-3 

3-4 

2 

jm 

0.897 

0.877 

0.897 

0.826 

0.856 

0.826 

0.856 

3 


.867 

.857 

.827 

.785 

.836 

.790 

.774 

4 


.857 

.777 

.797 


.829 

.703 

.741 

5 




.777 




.723 


Based on these R values, figure 3-3 (S = 2, r = 3; T = 4, [2 retained 1]) seems 
to be an optimal rotation pattern since it yields the minimum R (R = .777 if 
Y = 2.01, R = .703 if Y “ 3.38). This indicates that, with this rotation 
pattern, the cost of AgRISTARS could be reduced to about three-fourths. 

A 

3. 5. 4. 2 Values of R for Rotation Patterns with Return for Various Values of y 
I n a previous section, an optimal rotation pattern was chosen based on the 

A 

values of R confined in two specified values of y (2.01 and 3.38). Now the 
empirical optimization will take place based on values of R computed from 

A 

various values of y* Table 3-9 presents those values. 
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TABLE 3-9.- VALUES OF VAR a, IN DIFFERENT ROTATION PATTERNS WITH 
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3. 5.4, 3 The Optimal Rotation Pattern 

The comparison of the values of R in table 3-10 indicates that the rotation 
patterns in figures 3-3 and 3-4 are those better than figures 3-1 and 3-2. 


Between figures 3-3 and 3-4, if the values of y are greater than 1,78, the 
rotation pattern in figure 3-3, that is, (S « 2, r " 3, T = 4, [2 retained 1]), 
is the optimal rotation pattern. 

This comparison can be figured by using the illustrations shown on the follow- 
ing pages (figs, 3-19, 3-20, and 3-21) in which the minimum values of R shown 
in figure 11 are compared with those of figures 3-1, 3-2, and 3-4. 

Figure 3-19 shows the comparisons among four figures, figures 3-1 through 3-4, 
in terms of R values. For T = 2, figure 3-1 is better than the others. 

Figure 3-20 gives the comparisons among three figures, figures 3-2, 3-3,. and 
3-4, for T = 2, 3, and 4 years. In figure 3-20, the R values in figure 3-2 
are often higher than those of figures 3-3 and 3-4. So, figure 3-2 is 
eliminated to provide figure 3-21 which gives the comparisons between figure 
3-3 and 3-4. (In fact, figure 3-4 needs one more year, T = 5, in order to 
achieve its full cycle). For T = 4, except for the period around T = 3, most 
of R values of figure 3-4 are higher than those of figure 3-3. Moreover, 
figure 3-3 represents a full -cycle rotation pattern (S « 2, r = 3, T = 4, 

[2 retained 1]). 

The length of time (T = 4 years) to complete the full -cycle for the rotation 
is shorter than that of figure 3-4 (T = 5 years). Therefore, figure 3-3 is 
more favorable than figure 3-4. 

Henceforth, figure 3-3 represents the optimal rotation pattern for T = 4. 
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Figure 3-19.- R values in figures 3-1, 3-2, 3-3, and 3-4 by y and T = 2, 3. 
u°r ^ £ ?? ^"iQure 3-4 yields a smaller R value than do the other figures, 
but of all four figures, figure 3-1 shows the smallest R value. For T = 
figure 3-1 shows the smallest R value. 
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Figure 3-20.- R values in figures 3-2, 3-3, and 3-4 with respect to y and T 
for T = 4, figure 3-3 shows a smaller R value than do the other figures. 
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Figure 3-21.- R values in figures 3-3 and 3-4 with respect to y and T. 
or T » 4, figure 3-3 shows a smaller R value than do the other figures 



4. CONCLUSION 


4.1 TH E SAMPLE DESIGN FOR FCPF 

4.1.1 THE OPTIMAL ROIATION PATTERN FOR FCPF 

In the case concerning the applicability of the rotation sample design for the 
AgRISTARS FCPF project, this study Indicates that the rotation pattern (S = 2, 
r = 3; T = 4, [2 retained 1]) Is the optimal rotation pattern. 

It Is worthv/hlle to recapture this optimal rotation sample design. 

4.1.2 THE ROTATION SAMPLE DESIGN FOR FCPF 

Average wheat acreage per sampling unit Is the characteristic to be estimated 
by the rotation sample design for each ''current" year. 

The elementary units sampled are "p1)<els" of the earth's surface for many 
selected regions of some countries such as the United States, Canada, 

Cl.S.S.R., Australia, India, China, Argentina, and Brazil. 

The sampling units are "segments", which are clusters of pixels. The sampling 
unit size Is large enough that the variation of the wheat acreages of a 
particular segment from year to year is usually less than the variation of the 
wheat acreages of different segments within a particular year. 

Each wheat region Is stratified Into L strata which are homogeneous In terms 
of wheat density (i.e., wheat densities vary little within a stratum but may 
vary considerably from stratum to stratum). 

Suppose a stratum h, h = 1, 2, •••, L, Is comprised of Nj^ segments. Out of 
N|^, a predetermined number of segments, for example n^ = 1, 2, or 3, are 
chosen to be In the whole sample of study. 

The number n^, v/ill be the same for each year of study, but the segments will 
be chosen by the rotation pattern (S = 2, r = 3; T = 4, [2 retained 1]), as 
shown in figure 3-3 which Is repeated as follows. 
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Year number 


Segment 

number 



(S = 2, r = 3) Rotation Pattern [2 retained 1; T = 4]. 
The basic model Is 

\s “ “t * '’s ®ts 
for 

t = 1, 2, 3, 4 
T 5 4 

s = 1, 2, 3, 4 
where 


(241) 


= average true wheat acreage per segment In year t; at are fixed-year 
constants 

bg = true segment variables applicable to all years with the assumption 

1>S ~ «(“’ “b) 

e^s = composite-segment-error variable of segment s In year t, with the 
assumption ®ts ~ *^e) 

A^g = current year "direct" estimated (from satellite data) or "previous 
estimated" wheat acreage of segment s of stratum h In year t 


In matrix form, the basic model, which Is shown In equation (241), can be 
written as 

a = Xa + Ub + le (242) 

r>j fsj fNrf ' / 
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where 


10 0 0 
0 0 0 1 


0 1 


0 0 


1 0 
1 0 


( 243 ) 


i 0 0 0 
10 0 0 
0 10 0 
0 10 0 
0 0 10 
0 0 10 
0 0 0 1 
0 0 0 1 


a = (aj^, «2» «3> “4 ^ *5 “4 is to be estimated as 


b = (bp bg, bg, b^)' 0, la- 


e = (ej^p e^p ej^g* ®22’ ®23’ ®33’ ®34’ ®44^ ~ ^°e 

~ ■ ^41’ ^12’ ^22’ ^23’ ^33’ ^34’ ^44^' 


Then, the estimate is deduced from the vector 


a = (X'H"^X)"^X'H"^a 


which is the BLUE of a, with 


„ " 1 + 2y /i ^ 1 2 \ _2 

^ ^ \ 1 + Y 1 + Y/ ® 
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where 


_ *2 hi 

- %r^ 


as estimated from the pilot data 


i-1 = 1 


T'+Tf 


1 + Y 

-Y 

0 

0 

0 

0 

0 

0 


-Y 

1 + Y 
0 
0 
0 
0 
0 
0 


H = I + yUU' 
0 


(251) 

(252) 


0 

1 + Y 

“Y 

0 

0 

0 

0 


0 

0 

-Y 

1 + Y 
0 
0 
0 
0 


0 

0 

0 

0 

1 + Y 

' -Y 
0 
0 


0 

0 

0 

0 

-Y 

1 + Y 
0 
0 


0 

0 

0 

0 

0 

0 

1 + Y 

-Y 


0 

0 

0 

0 

0 

0 

“Y 

1 + Y 


(253) 


X'H“^X = Circ[2(l + y)» -Y, 0. -yI 


(254) 


Therefore, the estimate of the stratum h wheat acreage at the year T h 4 will 
be 


^4h ^h“4 


with the variance 


Var 




\ Var 


(255) 


(256) 


And, if is the current-year stratum-h average wheat yield estimate, the 

estimate of the total wheat production in the wheat region of interest will be 


Prod. = £ A, 


"4 " ph ^'4h; 


(257) 


From table 3-10, note that when y =* 4.00, R = .68. That is, with this optimal 
rotation sample design, the reduced variances of the rotation (multiyear) 
design estimates go to 68 percent of the corresponding variances for the one 
year estimates. 
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4.2 LIMITATIONS 


The first limitation, as mentioned previously, is that the target population 
must yield nj^ = 1, 2, or 3 sample segments per stratum h. This kind of target 
population was named Group I, 

The second restraint is that the sampling unit size needs to be large enough 
so that the variation of wheat acreages of a particular segment from year to 
year is less than the variation of wheat acreages of different segments within 
a particular year. 

Hartley et al. (ref. 4) have shown that if the estimate of a's (hence of y) is 
consistent, a formula similar to equation (29) applies asymptotically. There- 
fore, the thir^i requirement is that the estimate of y needs to be consistent 
so that the formula shown as equation (250) is valid. The estimate y is 
obtained from the pilot data. 
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APPENDIX A 


THE INVERSE OF A SPECIAL MATRIX 


For a special matrix of the form 


(a,j) . 


it is easy to prove that 


a b ••• b 
b a • • • b 

,b b ••• aj 


nxn 


la,j| = (a - b)""l [a + (n - 1)'’] 
and that the inverse is also of the same form, but now 

.ii a + (n - 2)b 


[a + (n - l)b] (a - b) 


and 


(A-1) 


(A-2) 


(A-3) 




-b 


Ca + (n - l)b] (a - b) 


i j 
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APPENDIX B 


INVERSION OF A SPECIAL TRIDIAGONAL MATRIX 

A special method which has been devised using finite difference calculus will 
be utilized to Inverse the following special tri diagonal matrix. This matrix 
Is of the form 

r— 

a c 0 0 • ♦ • 

c b c 0 

< 

0 c b c 


(B-1) 


c b c 0 
0 c a c 

The last column of the Inverse of matrix (B-1) will be denoted as 

V 


X, 


The equation of the definition of the matrix Inverse 

AA“^ = I 


(B-2) 


will give a system of T equations. Not considering the first and the last 
equations, the (T - 2) middle equations are 

^1-1*^ + x^b + x^^j^c = 0 (B-3) 

which are 

X 1.1 - 2x, f (B-4) 

denoted as 

a"x. = xx^ (B-5) 


B-1 


. C3 ftr|k 


where 


A solution Is 


A 5 -(b + 2c)/c 




so, ecjuation (B-4) becomes 


- (2 + X)p + 1 a 0 

which has two roots p^ and pg. 

Hence, a general solution will be 

» ap\ + 3pJ 


With che condition 


4X + > 0 


(B-6) 


(B-7) 



(B-8) 

t equations which are 


BpI) ■ 0 

(B-9) 

+ spj) - 1 

(B-10) 


(B-11) 


which Is the positive discriminant of equation (B-7), the two roots p, and p, 
wi 1 1 be i 2 


Pl» P2 - 1 + Y * 

Two equations (B-9) and (B-10), hence, will give 


3 =* 
a = 53 


^ 4 '')] * - pS]|-' 


(B-12) 


(B-13) 


B-2 


where 


Therefore , 


3 -(cp2 + ap2)/(cp2 + apj) 


*T ■ “"I * 


(B-14) 


(B-15) 


B-3 
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APPENDIX C 

THE INVERSE OF A CIRCULANT MATRIX 

The following presentation Is provided by Dr» H, G. Newton, professor, 
Texas A«M University. 


A clrculant matrix Is a matrix of the form 


A n 


which Is denoted as 
Equation (C-l) can be written as 
where 

P =* (Pp — » Pt^ 

A =• dlag (Xj, Ay) 
w1 th 


”a(l) 

a(2) 

a(T) 

a(T) 

a(l) 

a(T - 1) 

_a(2) 

a(3) 

a(l) 

» ClrcCadl, 

•••» a(T)3 


A » PAP' 


X 12Tr(j-l){k-l) 

A, 5 £ a(k)e T ; .1 = 1 , 2 , •••, T 

k*l 


T^/^ P. = 


1 


-12ir(j-l)a)/T 

- 127 r(j-l)( 2 )/T 


- 127 r(j-l){T-l)/T 


(C-l) 


(C-2) 


(0-3) 


; 0 = 1, 2, •••, T 


(C-4) 



since 


P'P « PP' = L 


(P')"^ = P 

P"1 = p' 


thus, 


A’^ = (> A P')"^ 

= (P‘ )"^ a“^P"^ 


= P a"^ P 


And, the (j, k) element of A"^ will be 


(A"^)jk ” ^ times k^*^ column of (a"^P') 


= .i; Pji(A-^P') 


1=1 

T 

E 

1=1 

T 


Ik 


= J2 a"^ times k**^ row of P) 


“ ^^I'^kl 


_ 1 V 6 
~ T 2 ^ 

' 1=1 


T -i2Tr[(j-l)(l-l) + (k-l)(l-l)]/T 


1 


I ^ g-i2Tr(l-l)(j-l + k-l)/T 
= 1 


T ^ 

' 1=1 


(C-5) 


(C-6) 


C-2 


hence 


and 


For example, 


, 1 i ,-12.(1-l)2(J-l)/T 
" T ' 


1=1 


(A"^) 


_ 1 ^ 1 
^ 1^ 


1 1 X ^-i2ti{1-1)2(T-1)/T 

('' >tt = T ,i: ^ 


A = C1rc(2(l + y)» -y» 0» -y) 

X, = 2(1 +y) - Y‘**0”Y~2 

JL 

Xg = 2(1 + y) + (“Y) + (-Y)e"'*®’^'^^ 

= ,2(1 + y) , • 

X3 = 2(1 + y) + {-Y)e“''^’''''^ + 

= 2 + 4y 


X^ = 2(1 + y) + (“Y) e 
= 2 + 2y 


-i6iT/4 


+ (-y) e 


-il8Tr/4 


,.-l, _l(l , I 

>44 -T 


1 ”2 
1(1 1 


3 ^4 ) 

1 1 I 


4 ( 2 2(1 + y) 2 + 4y 2 + 2y j 


-1 1 + 

• ■ 8 r ‘ 1 + Y 1 + 2 y 


(C-7) 


(C-8) 

(C-9) 


(C-10) 


(C-11) 


C-3 


